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DESCRIP TION OF THE INVENTION 
Definitions 



Acoustic, radiation prediction : is an engineering process that aims at 
predicting the sound pressure level in the surrounding field of a vibratinq 
structure (see figure 1). Surface vibrations create acoustic pressure waves 
that propagate to the receiver. 

Inverse numerical acoustics : is an engineering process that aims at derivinq 
the surface vibration velocity from field sound pressure measurements (see 
figure 2). As such, it is the inverse problem than acoustic radiation prediction. 

Discretisation Methods .- in order to solve acoustic radiation prediction 
problems and inverse numerical acoustics problems, computer models are 
built using numerical methods. The most frequently used numerical methods 

mfth^H .™/? c a a!,f ,e -S ent meth0d (BEM) and finite-infinite element 
method (FEM/IFEM). These methods are relying upon a discretization of the 

2? main (FEM/IFEM > or of its boundary (BEM). and are solving 
numencally the wave equation, in the frequency domain. There exist 
substantial published literature concerning both theoretical and practical 
aspects of these numerical methods. The major computational effort involved 
within the application of such numerical methods concerns the solution of a 
large system of linear equations. 

Background of »h q invention 

The acoustic performance of manufactured products is becominq an 
extremely important aspect in the design and development process, not only 

?n im ^%?l?? m f° rt ? the USer ° f the product < e -9- P^ngers in a car, in 
an areraft). but also to reduce the nuisance to the surroundings (e g 
habitations close to highways, to an airport). { 9 " 

The invention is related to two computer-aided engineering processes that 

the ab^ r tn Va,U H w 9 + and ° ptim ' 2ing the acoustic Performance oSSSiSE 
l»t 8 5? ° prBd,Ct the acoustic radiation pattern of a vibrating structure 

P^ictioT a~d m ?h? d a hnLT aSUred SUrfa S e Vlbrations < acous « c radi ^« 
2^S™ t P to r9cover surface vibrations onto a vibrating 

structure from measured field sound pressure level. The latter one is 

t °V y manner tQ perfbrm 80Ure8 identification, when it is 
rolling tire) * measurement devlces onto *e structure surface (e.g 

wwrh th!! 9 ™^".? processes are relying on numerical methods, amongst 
2£E£J?t Th I P^ PU L ar 3re tne various fo rms of the boundary 

element method (the direct mono-domain, the direct multi-domain and the 

S^uKES? ^ e,em8nt method ' that 030 be tended to 

handle unbounded regions, e.g. using infinite elements. 

The computational cost associated with numerical methods for acoustic 
performance prediction is usually very large, and is linearty proportional to: 
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• the number of operational conditions (may be about a 100 different 
cases); 

• the number of frequency lines to be evaluated for obtaining a 
representative response function (will typically be about 100 to 200); 

• the number of design variants to be studied. 



Acoustic Transfer Vectors 

The invention relies on any of those numerical methods, and is based on the 
known concept of Acoustic Transfer Vectors (ATV's). ATVs are input-output 
relations between the normal structural velocity of the vibrating surface and 
the sound pressure level at a specific field point (see figure 3).. ATV's can be 
interpreted as an ensemble of Acoustic Transfer Functions from the surface 
nodes to a single field point. Literature refers sometimes to that concept as 
contnbution Vectors or Acoustic Sensitivities. 

A 7\!f on, y de P end on tne configuration of the acoustic domain, i.e. the shape 
of the vibrating body and the fluid properties controlling the sound 
propagation (speed of sound and density), the acoustic surface treatment the 
frequency, and the field point. They do not depend on the loading condition 

I.K.i^ ept ATS f S ' and their P ro P ertjes . is not "ew, and has already been 
published in several scientific papers. 



Summary of th» invention 

The present invention will mainly be described with reference to acous'-- 
waves but the invention is not limited thereto. The concepts underlyinq the 
present invention may be applied to any vibrating energy wave which can be 
1 J!? 1 ?*? by , W3Ve ec > uations - For example, the present invention may be 
2L«!!3 ? e,ectn?ma 9 netic waves. The present invention may provide v 
method for computing a vibrating wave Transfer Vector based on the 
reciprocity pnnciple, comprising the steps of: 

- simulating positioning of a monopole, omnidirectional wave eneray source 
at a reference position remote from a body; 

- computing a boundary vibration amplitude of the wave generated bv the 
source at a surface of the body; y 

Vector 9 fr ° m the boun *»'y vibration amplitude said vibrational Transfer 

?nSnl^ v S 87 a ' S ° indUde 3 method for ^Puting an additional 
viorational Transfer Vector compnsing the steps of: 

computing at least a first and a second vibrational transfer vector at a first and 
a second predetermined frequency, respectively, and 
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computing the additional Transfer vector at a frequency intermediate the first 
and second frequency by interpolation between the first and second tranfer 
vectors. 

The present invention may also include a method to compute a Modal 
Acoustic Transfer Vector (MATV) from an acoustic transfer vector (ATV) for a 
modal space defined by the Eigen-frequencies of a body. comDrisina the 
steps of: f v 

projecting the ATV into the modal space, and 

using the MATV to predict a response of the body or the effect of such a 
response at a reference point remote from the body. 

The present invention may also include a processing engine adapted to carry 
out any of the methods of the present invenrtion. 

The present invention also includes a computer program product for 
executing on a computer, the computer program product executing any of the 
methods of the present invention. The computer product may be stored on a 
computer readable data carrier such as one or more diskettes or one or more 
CD-ROMS or may be transmitted voia a telecommunications netv/ork such as 
the internet 

The present ibvention also includes a method of inputting at a near terminal a 
representation of a body and coordinates of a reference point and transmitting 
these to a remote terminal running a program for executing any of the 
methods of the present invention, and receiving at a near location an output 
of any of the methods. The output may be one of: 

an ATV. a vibrational amplitude such as an acoustic pressure level, a surface 
vibration of the body, a revised design of at least a part of the Lx)dy. 

The invention consists in one aspect in a new methodology to dramatically 
speed up the evaluation of the Acoustic Transfer Vectors and their 
exploitation. In particular the use of the reciprocity principle allows to improve 
the accuracy of the ATV's because they can be computed and stored on an 
SS 6 ? ^asis, and to reduce *e computational effort, since the solution 
needs to be performed for a number of excitation vectors equal to the number 
of field points of interest, instead of to the number of boundary nodes. 

Further on, the use of efficient interpolation techniques is very appropriate for 

SJtnw tK 9 - ?' dUe to their smootn frequency dependence, and allow to 
restnct their evaluation to a restricted set of master frequency lines, which are 
then de-coupled from the excitation frequency contents. This results in 
extreme gains in computing time and data storage space. 

The developed methodology results therefore in a reduction of the 
computational cost attached to acoustic performance prediction, makinq it 
quasi insensitive to the number of operational conditions, to the number of 
frequency lines of interest, and to the number of design variants. As such it 
allows the acoustic performance prediction tool to be driven by an automated 
design exploration and optimization tool 
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Benefit of the invention 



Acoustic performance prediction is applied to many different kinds of 
products, in a process that is well standardized. As an example, let us 
So^tf' 3n a " tomotive Powertrain structure, whose goal is to minimize the 
acoustic sound emission, and try_to_quantify the benefit due to the invention. 

A reasonable boundary element model will include about 7000 
nodes/elements. It should be analyzed under 50 different operating 
conditions, corresponding to different engine speed (RPM's). For buildina a 
ftjH frequency response, about 200 frequency lines need to be computed 
The performance at 10 microphone locations is of interest, for example to 

procedure) radiat6d S ° Und P ° Wer (aCCOrX,in9 to tne ^MtJmIb? 

'^SkSSSfT^l ^ f ° r P,? cessin 9 such * "i°del on an engineering 
fi££f££Z, aSSd . ,° n followm 9 fi 9 ures: the matrix assembly and 
??2F? 0n PmCeSS 18,588 5 ' 400 secon ds per frequency, and the back- 
substitution process (needed for each excitation case) takes 3 seconds oer 

IS3S!?- ^ S H min9 M 3bOUt 20 diff6rent desi 9 n aitemaLL n^ed to P be 
analyzed in the des.gn exploration process, we get following estimations: 

* ^drtional acoustic prediction technique, we can estimate the 
£? timf 1 ? ° nal 1 8ffl ? rt t0 56 the number of frequency lines (200) times 
SsesT VZ?tT 9le *f*™!* ana ' v S* (5.400 sec 7 3 sec \ 50 lold 

tZZSL cTrrUn* H Umb , e o£ d6Si9n . altenia « v es (20). |.e. 22.2 million 
seconds, corresponding to 257 computing days. 

• Using the Acoustic Transfer Vectors technique, computed in a direct wav 
the model has to be solved for a number of excitation cases equll to The 
number of nodes (7000). but since the exploitation time (scalar^rodi^t il 

SSfSLS^ ™ atf ° nal 0051 '^Pendent onlhe Ser of 
design alternatives, and of operating conditions. Therefore the total 

tZsteT^L^ iS es f imated to be the of frequenc^iines (20^ 

times the time for a single frequency analysis (5,400 sec + 3 sec x 7000 

vectore) 528 mi,,ion seconds ' ^ S poVd^ X to 7 °6? 

• Using the invention, i.e. the Acoustic Transfer Vectors techniaue 

™ZTe ZbTo ™T? 00mblned With « wSoTbS" 

in « 1 ^ Ster fre q uen *es, and neglecting the time for ATV 

2fartSS2 n th? ,d e ^ p,oit ? tion - W V est »' mat * the total computational 
effort to be the number of master frequency lines (10) times the time for a 
single frequency analysis (5.400 sec + 3 sec x 10 load cases) | e 5^ 300 
seconds, corresponding to 0.6 computing days. 
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THEORETICAL BACKGROUND 



Notations 



A, B t H. Q Fluid domain matrix operators 

A TV Acoustic Transfer Vector 

ATM Acoustic Transfer Matrix 

c Sound velocity 

C Coupling matrix 

G Green's function 

/ Frequency 

* Wave number 

K. D, M Fluid stiffness, damping and mass matrix operators 

J Imaginary number 

* Surface normal 

P Acoustic pressure 

P*> Boundary acoustic pressure 

S Surface of the vibrating (emitting) structure 

u - v Complex matrices (SVD) 

v Acoustic particle velocity 

v * Boundary structural velocity 

v * Structural velocity 

x Observation point 

y Source point 

0 acoustic admittance (complex) 

P mass density (complex) 

<* Single layer potential; singular values 

H Double layer potential 
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<o Circular frequency 

Internal degrees of freedom (superscript) 
J Interface degrees of freedom (superscript) 

Governing Fqnatirtrf 

Acoustic wave propagation in bounded or unbounded regions is qoverned bv 
the wave equation. Assuming an harmonic behavior, the time and space 
vanables may be separated, leading to Helmholtz' equation: 

V 2 p + k 2 p = 0 
where k is the wave number ©/c. 

Four numerical formulations are available for solving Helmholtz equation over 
k ^?!? n U ? n S. ed by a surface * with P ro P*r boundary condftons and to 
tUe1n^on° ranSfer VeCt ° re ^ imPr ° Ved Way ^nciptes of 

1. The direct boundary element method (DBEM) formulation usina the 

Tn be us P ed r^^hT^ " "T^* Vanab ' eS ' This ^^ion 
can De used to model homogeneous intenor or exterior domains. 

2 " I™ 1 m . u,ti - domaln dir ect boundary element method (MDDBEM) 

SSS?° n, TH U - 8,, } B thS 3COUStic pressure and ve,ocitv as boundary 
vanables. This formulation can be used to model a non-homoaeneous 

JST l y diV,d : ng * in 3 set * "omogeneous sub-domains thaTarL Sed 
together by continuity conditions. 

3 * sinoli"^^ n( ? ary e,ement method (,BEM > formulation, using the 
single and double layer potentials (acoustic velocity jump and nressurp 
ump respectively) as boundary variables. This forrUaton can be used 

SdSte .n ho i mo 9 1 eneous do ^in. where both sides of the boCndarJ 
radiate. In particular, it can be used to model an interior domain an 
extenor domain or both simultaneously. °omain, an 

4. The acoustic finite/infinite element method (FEM) formulation usina the 
nS!?2 PreSS K re 38 b ° Undary Variables - This formulation S'be* sld to 
SLnel^y 0 ^^ interi ° r d0main ' a " <~ *>™in - tth 

Direct Boundar y Eleme nt Method (DBErV» Farmnlati^n 
Using the DBEM formulation, the acoustic pressure at anv noint nf 
homogeneous fluid domain containing no acouSic source can ^beexSresLd 
deSefllf 6 3COUStiC PreSSUre °" me b0Undaf V domamand 6 
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(1) 



3 y 



where p(?) is the acoustic pressure on the boundary and its normal 

dn y 

derivative, n y is the inward normal at point p on the boundary and G{x\y) is 

the Green's function. Making use of the Euler equation, equation (1) 
becomes: v ' 

Pi*) = jp(y) ^^ JSy + JP«> Jv(y) G(x\y)dS y (2) 



where v(y) is the normal acoustic velocity on the boundary, p is the fluid 
mass density and a, is the angular frequency. Note that the boundary 
rSat!onshTp IOCity ' S re ' atGd to the structural velocity through the following 

= v, GO + P(y)p(y) (3) 

where v s iy) is the structural velocity and p(y) the boundary admittance. 
Using equation (3), equation (2) becomes 

fcW^^ dS y + >P«» fmp(?)G(x[y)dS y +7 po> f Vs (y)G(x\y)dS y 



(4) 



ISu^h 1 tn if m ? e domain and on rts bou ndary. Nevertheless, when 
evaluated on the boundary, the Green's function and its normal derivative 
become s.ngular. Whereas the last two integrals of equation (4) are regular 
vaTue^ns? ,S S ' ngU,ar and shou,d be evalu ated in the Cauchy's principal 

c&P&-PV.jp(?)?^ d S y +JPG > k(y)p(j,)G(x\s)dS y L(y)G(x\y)dS y 
s y s y 

(5) 

where 
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in the three dimensional space. Note that c(x) = i/2 for a smooth surface 
around x. 

Once discretized using boundary elements and evaluated at the mesh nodes 
equation (5) leads to the following matrix system: 

Mu^l-MM (7 ) 
where the subscript b stands for boundary. Similarly, equation (4) gives: 

P = {<*} T {p b } + {>nVM (8) 
Combining equations (7) and (8) leads to: 

p = {atv} r {v b } ( 9 ) 
where {atv} is the Acoustic Transfer Vector, given by. 

w T -w T [^r'w+w T do 

The Acoustic Transfer Vector (ATV) is therefore an array of transfer functions 
between the surface normal velocity and the pressure at the field point 
Finally, when the pressure is evaluated at several locations, equation (9) can 
be rewritten as: v ' 

{p) = [ATMYM (U) 

where the Acoustic Transfer Matrix [ATM] is formed by the different Acoustic 
Transfer Vectors. 

. Multi-Domain Direct Bo i mr f arv Element Method (MDDBFM) Formulation 

In case the acoustic region includes partitions characterized with different 
fluid matena properties, the above formulation can not be directly used since 
it is only valid for homogeneous domains. In such cases, a multi-domain 
integral formulation has to be used. The global acoustic region is then 
decomposed into sub-domains, with the requirement that within a sub- 
domain, the fluid properties need to be homogeneous. At the interface 
□etween the sub-domains, continuity conditions are enforced. 

For sake of simplicity, let's consider here a two-domain model. Equation (7) 
for the first sub-domain can be written as: 



4" Af 

L Af A?\ 



fcJH? 3M 
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11 



where superscript / stands for internal degrees of freedom and superscript j 
stands for interface degrees of freedom. Similarly, equation (7) can be 
rewritten for the second sub-domain as: 



(13) 



The continuity of the normal velocity and acoustic pressure has to be satisfied 
at the interface between of the two sub-domains: 





4 


K 








At 








At 


B*_ 




N =-V 2 = v 

{pi = p{ = p 



(14) 



Combining equations (12), (13) and (14) leads to the global system of 
equations: 3 



4" 






A Ji 






0 




*i 


0 


4 


B» 



0 
0 



Pi 



P'z 



*t 
0 

0 



0 
0 
B$ 

*tl 



ft} 



(15) 



U . S k 9 , eq " ation ( 15 )' one can define the Acoustic Transfer Vectors for the 
global system: 



= {atv] 



ft} 



(16) 



and similarly derive the Acoustic Transfer Matrix [atm] when the acoustic 
pressure has to be evaluated at different locations. 



Indirect Boundary Foment Method (IB EM) Form.. 

rin£L a ,H in K lreCt formu J ation ' acoustic pressure at any point of the 
domain can be expressed in terms of single and double layer potentials [2]: 



<• y ^ 



(17) 



where o{y) and My) are the single and double layer potentials, respectively 
k w i! euman P rob,em Ce. structural velocity boundary condition) and 
when both sides of the boundary have the same velocity (i e SSn sheH this 
equation can be expressed only in terms of double layer potentials 
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s y 

Equation (18) can be evaluated on an arbitrary surface S' . It can be derived 
with respect to the normal at xeS\ pre-multiplied by an admissible test 
function fy(x) and integrated over the surface S' : 

Now, let the surface S' tend to Sand make use of the Euler equation: 

- jja^v(x)SMx)dS = ^My f^^ &(Z)dS y dS x (20) 



S S S 



After discretization using boundary elements equation (18) becomes: 

P = {d} T \M (21) 

Similarly, equation (20) becomes: 

feW = Wv 4 } (22) 
Combining equations (21 ) and (22), we obtain: 

P = {atv} T {v b } (23) 
where {atv} is the Acoustic Transfer Vector, given by: 

w T -w T ferw (24) 

When the acoustic pressure is evaluated at several locations, equation (23) 
can be rewritten as: 

W = W(v 6 } (25) 

where the Acoustic Transfer Matrix [ATM] is formed by the different Acoustic 
Transfer Vectors. 



Finite/infi nite element method fFEM/IFEM) formulation 

The solution of Helmholtz" equation based on a finite/infinite element 
approach leads to a system of linear equations: 
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t^I+^Z)]-^ [m^p] = -jpo^F) (26) 

where K is the fluid stiffness matrix, D is the fluid damping matrix, M is the 
fluid mass matrix and F is the forcing vector, defined by: 

(27) 

where C is a coupling matrix. 

Combining equations (26) and (27), ), we obtain: 

{p) = [ATMYM (28 ) 
where the Acoustic Transfer Matrix [ATM J is given by: 

WJT^pP^JmI (29) 

Finally, for each finite element grid point, the acoustic pressure can be written 

P=>{fUv} T M (30) 

where {atv} is the Acoustic Transfer Vector, given by the corresponding row 
of the (ATMf matrix. a 

Use of the Recipr ocity Prin^i plo 

^r ke ^ m « nt in ihe described invention is the use of the reciprocity principle 
for efficiently computing the Acoustic Transfer Vectors As shown n 
equations (10£ (15), (24) and (29). the evaluation of any field T^oi^elated 
ATV requires the tnveraon of a matrix, whose order is approximately eaual to 
the number of nodes of the discretized geometry (mesh) V Q 

Jo h |vin n a « e rj?tL eVa l U f ti0n ° f the inverSe of an order n is performed by 

tnZVL^V? l,near e 9 uatlons . whose right-hand-side is built from the 

S&S^Jg*" 9 ***** t0 ° ne mat ™ ^orization. «d n back* 

[AIX]=[I] (3l) 
For large matrix orders, the computational effort and time reomrpri for 

sat^^ «*«» - «w sm«*«$?f a ^z? n £ 

Looking to equation (10). we can deduct that an ATV is a vector for which 
each coordinate corresponds to the pressure at the conlslo^g Z^SnX 

3 h" exc,tat,on (vibrating velocity) at one poirKlSSi^^E 
surface, and no excitation elsewhere (see figure 4). '™ eie menr on the 



03-07-00 11:29 32 16 480528 P. 17 R-592 Job-GI6 

03-07-2008 12:55 SiRD & GOEN & CO 32 16 480528 P. i 7/35 

14 



arv,. = {a/vH 1 1 (32 ) 



According to the reciprocity principle, the pressure response at a first location 
due to an excitation at a second location is strictly equal to the pressure 
response at the second location due to the same excitation at the first 
location. Therefore, if a source of strength Q, at position 1 causes a pressure 
p 2 at position 2. and if a source of strength Q 2 at position 2 causes a pressure 
Pi at position 1 , we have the following relationship: 



or = 



(33) 



Using the reciprocity principle to evaluate the ATV related to a specific field 
point means that a monopole source (point source) of strength Q, should be 
positioned at the field pcint location 1 , and the pressure is calculated or. th-=» 
boundary p 2 (see figure 5). 

Equation 11 explicates how the ATV will be exploited to convert the boundary 
surface velocity (representing vibrations) into the field point acoustic pressure. 
Since a structure may be complex, containing T-junctions and other 
discontinuities, surface vibrations are better defined on surfaces, i.e. on finite 
elements, rather than at grid points (where the normal directly is not always 
uniquely defined). Therefore, it is required to define element-based acoustic 
transfer vectors. 

In the case where the source placed on the boundary (in 2) is distributed over 
a small surface, equation (33) may be rewritten as: 

PiQx = [p*v 2 dS (34) 

s 

U^S? V? the source ve,ocit y- In the particular case where the source is 
distnbuted over a single element, the pressure at any point on the surface can 
be approximated as: 

Pz^iPiW (35) 

where {N} is the interpolation function and <p 2 ) is the pressure evaluated at 
tne element nodes. Similarly the source velocity can be approximated as: 

v 2 =(N){v 2 } (36) 

where fvi} is the source velocity at the element nodes. Equations (34) to (36) 
lead to the following relation: ' 
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P&^faftcM (37) 
where fCJ is the element coupling matrix defined as: 

[C e ]= [{v}(N)dS (38) 
For a monopole source, the source strength is given by: 

V\ — A (39) 
where A is the monopole source amplitude. Combining (37) and (39): 

which can be rewritten as: 

Pi = [ ATV ,\v 2 ) (41) 
where [ATVJ is the element-based acoustic transfer vector, computed by: 

Applying at the field point a monopole source with a unit amplitude (A=J and 
©-2^ we finally obtain the following expression for evaluating the element- 
based acoustic transfer vector 

[ATV € ] = M-( Pl )[ Ce ] (43) 

showing that the ATV related to a given field point can then be computed from 
tne Boundary sound pressure values, correctly scaled. 

The use of the reciprocity principle for evaluating the ATV's leads to two 
fundamental advantages: 

1 " Imt^T 5 0 ?" bS * d iL eCtIy com P uted on an element basis, removing the 
ambiguity inherent to nodal values, in case of sharp edges and T- 
junctions, where the normal vector is not uniquely defined. 

2 I^f^ re, ^ ed to a given fie,d point 030 then be computed with a single 
excitation vector, representing the monopole source excitation at the field 
point. Thjs leads to a single back substitution step. Since the number of 
field points (microphones locations) is usually several orders of magnitude 

o a 7n S r i^o m n e * n " mber ° f boundary nodes, this results in tremendous 
gains in computational effort. 
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Interpolation techniq ue 

i" ° f S ° Und Wave P ra P a 9a«on in an open space, the fluid domain does 
exhibit any resonant behavior. The sound field is fully determined by various 

Zf^T^T^ 5 ' SU ° h " "f eCtion by soIid ob j ects - absorption on damping 
r **ZL - w^? Ct, ° n amund and wave interference. Therefore. Acoustic 
™2S Vectors are rather smooth functions of the frequency, and ATV 
c£?~ J*T 030 be accurate, y evaluated at any intermediate frequency, called 
slave frequency, us.ng a mathematical interpolation scheme, based on a 

nn^fZTJ r ? fre , qu ! n deS ' 03,1611 master frequencies. It is important to 
am £££ «™ ct " ra, . v ' brat, ° ns be similarly interpolated, since these 

are directly dependant on the highly resonant dynamic behavior of any 
structure (vibration mode shapes). ur any 

Two interpolation mechanisms are used: 

1. A linear interpolation scheme (see figure 6); 

2. A cubic spline interpolation scheme (see figure 7), which is particularly 
well suited for approximating ATV's. ^"'cuiany 

Other interpolation schemes can be used. 

JndaTenS SCHeme ^ GVa,Uatin9 the ATV ' S ,eads to fo "™"9 

1 . It allows to build a complete frequency dependant ATV. with a very limited 
ZnuS ° f frequencv to be computed. In order to build a comp ete 
respons f ™ th a traditional approach, the number of frequent 
lines to be computed is usually over a hundred. Using an interpolation 

° n V ? W maStef frequen ^ " nes " eed * * b^TpSy 
computed, any other frequency being interpolated. The gain in computing 
time is usually over an order of magnitude. computing 

2. Similarly, the ATVs have to be stored in a file database in order to be 
further exploited. The interpolation mechanism allows to store only t£ 

EK^SET* ^ in th6 d3tabase The "^spending verfor at a™ 
mterniediate frequency may then be directly evaluated at a neglectible 
cos and therefore does not need to be stored in the database The ga n 
in disk space is usually over an order of magnitude. 9 

3. For some specific applications, such as automotive powertrains the 
acoustic performance has to be evaluated for many different operational 
condraons (e.g. different rpm's). each of these having a^S^ S^e^S 
contents. Using the ATV interpolation mechanism, the same set of mSS^ 
frequencies can be used to build a unique basis of ATV^hfch are thSi 

^ d a f a c ° mbi "n With the tAaunl -spons^^sp'ondtnrto 
op^otr^ndi2on any ' ntemied ' ate *r th P e specific 



03-07-00 11:29 32 16 480528 P 20 R-592 Job-616 

03-07-2000 12 = 57 BIRD & GDEN & CO 32 16 480528 P. 20/35 

^ 17 



Inverse Numerical Acoustics 

Equations (11) and (25) are the basic relations for the inverse numerical 
acoustics theory. The sound pressure is measured at a large number of 
microphone (fiefd points) and the boundary velocity can be obtained using the 
relationship: 

{v b } = [ATMY ! {p} (44) 

However, the inversion of the acoustic transfer matrix is not obvious, because 
the matrix is generally not square but rectangular and because it involves a 
Fredholm equation (of the first kind for the indirect approach and of the 
second kind for the direct approach) and is ill-conditioned. Therefore the 
matrix inversion is performed using the singular value decomposition (SVD) 
technique, which allows to sofve singular or close to singular systems. It is 
based on the fact that any n x m complex matrix can be written as: 

[ATM] = [v\oluY (45) 

where the superscript H stands for 'transpose complex conjugate 1 , [cr] is a 
diagonal mux(/i,m)x min(n y m) real matrix, [u] is a mxminfara) complex 
matrix and \r] is a n x min(n, m) complex matrix. The coefficients of [cr] 9 
called singular values, are stored in an decreasing order and matrices [u] 
and JV] are such that: 

[vnvMuTPMl] (46) 
From equations (45) and (46) it follows: 

[ATM]- 1 =[ul<xF[vr (47) 

Nevertheless, singular or ill-conditioned matrices contain null or small singular 
values. This results in infinite or vsry large values in [cr]~ J . Whereas the 
infinite terms clearly lead to an infinite solution, the very large values lead to a 
solution of equation (44) that will be very sensitive to the right hand side 
variations- Therefore, small errors on the RHS will result in very large errors 
in the solution. To avoid this problem, the small or null terms of [cr] are set to 

zero in [cr]" y . This process is called truncated singular value decomposition. 
The singular values are set to zero as soon as cr, < ctcrj where a is a 
tolerance parameter. 

Panel Contribution Analysis 

Panel Contribution Analysis (see figure 8) is a fundamental engineering tool 
to guide the product refinement in the development process. Equations (9) 
(16), (23) or (30) can be rewritten as: 
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nt _ 



(48) 

where /«> is the total number of elements. Subsequently, the contribution of a 
er e menL 9,Ven * SUmmation of e ' e ™* cSntnbu'uons "ve? ie panel 



where pe stands for panel elements 



(49) 



Modal A TV 

The engineering process to compute the structural velocity f Vh i on a vibratinc 
surface rel.es usually on the structural finite element meti?od and ften on 
modal superpo S ,t,on approach, where the structural response is expressed as 
a Imear combination of the mode shapes of the body: expressed as 

W=KW) (so) 

Combining (50) with Equations (9). (16). (23) or (30) leads to: 

P = {at V } T [o n \mpf} (51) 

where 

M r [o.Mm«*vf (52) 
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Optimization ^p^Ty 

nriU e I h0d h ? ? < ? n devel °P ed al, ° ws numerical acoustic radiation 
pred.ct.on calculation to be extremely fast, the main computational intensive 

T a hL ot CU, f t,0n ° f H 16 ATVs > beln 9 done as a Pre-processing step 
Therefore, it is directly possible and practical to integrate the acoustic 
radiation pred.ct.on within an optimization loop, where the objective function is 
the acoustic performance (sound pressure levels, or radiated Dower 
according to the IS03744-1981 procedure), and where the design variabtes 
are either structural or acoustic design variables. 

1^^^°" Pr ?°t SS f Shown at fi9ure 9 " 11 wi ' 1 tyPical'y combine a 
ma™ dynamic finite element solver together with an ATV-based (or 
MATV) acoustic prediction tool. 1 

Im plementation, 

I^n meth °? 3 ° f thG ^ Sent invention m ay be implemented on a processing 
eng.ne such as a workstation or a personal computer. The processing engine 
^wan? ase [ v f r accessible via a telecommunications network such as a LAN 
Sf^ ^1 Internet, an Intranet. The server may be adapted to carry out any 
of the methods of the present invention, for example a descriptor file of a 
body may be entered at a near terminal an transmitted to the server via the 
internet The server then carries out one of the methods of the present 

*r n if° n „f "S^ 3 near locatlon e "9- an e - mail box . any response of 
any of the methods jn accordance with the present invention Such a 
response may be, for instance, an ATV. a vibrational amplitude such as an 

tSSTpS^SS body! 8 V ' brati ° n ° f ^ * ^ ^» of * 

The methods of the present invention may be implemented as computer 

RnS m Th Wh ' Ch may be 8tored on data su ch as diskettes 2? a> 

ROMS These programs may also be downloaded via the Internet or anv 
other telecommunications network. y 
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CLAIMS 



1 A method for computing a Wave Transfer Vector based on the reciprocity 
principle; comprising the steps of: 

simulating positioning of a monopole, omnidirectional wave energy source 
at a reference position remote from a body; 

computing a boundary oscillation amplitude of the wave generated bv the 
source at a surface of the body; 

deriving from the boundary oscillation amplitude said Wave Transfer 
Vector. 

2. The method of claim 1 wherein the computing step is carried out bv a 
numencal method. 3 

3. The method according to claim 2 wherein the numerical method is the 
finite element method. 

4 * l^SS^^f^ ? Wherein the numerical m ethod is a combination of 
the finite and infinite element methods. 

5 " 7 ieth0 , d of daim 2 wherein the numerical method is the direct 
boundary element method. uireoi 

6 - £ma^ " UmeriCa ' meth ° d is ^ mu*. 

? ' Su^dT^ the nUmerical ~ - the indirect 

8 " InusU^u^ 0 ^ 9 *° fmkiUa Ci *' lm ' Wherein Wave source is a " 
9 the m steps d o? r C ° mput,ng an addi *°n*< Wave Transfer Vector comprising 

r S eco^n^t St 3 3 second wave transfer vector a t a first and 

a second predetermined frequency, respectively, and 

computing the additional Wave Transfer Vector at a freauencv 

Sr^S^J?* 3 ? S6COnd freqUSnCy by '"terpola^n beSn ^ 
first and second Wave Transfer Vectors. 

"'XtSXS?!^ *• in,en "" a,ion technk,ue ls 4,16 «* 
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12. A method to compute a Modal Acoustic Transfer Vector (MATV) from an 
acoustic transfer vector (ATV) for a modal space defined by the eigen 
frequencies of a body, comprising the steps of: 

projecting the ATV into the modal space, and 

using the MATV to predict a response of the body or the effect of such a 
response at a reference point remote from the body. 

13. The combination of methods of claims 1 and 9 for supporting the acoustic 
radiation prediction engineering process, in particular but not limited to 
automotive engine and powertrain noise radiation application. 

14. The combination of methods of claims 1, 9 and 12 for supporting the 
modal-based acoustic radiation prediction engineering process in 
particular but not limited to automotive engine and powertrain noise 
radiation application. 

15. The combination of methods of claims 1 and 9 for supporting the inverse 
numencal acoustics engineering process, in particular but not limited to 
automotive engine and powertrain noise radiation application. 

16. The combination of methods of claims 1. 9 and 12 for supporting the 
modal-based inverse numerical acoustics engineering process in 
particular but not limited to automotive engine and powertrain noise 
radiation application. 

17. The method of claim 12. wherein shapes combined to build a structurall 
response are Ritz Vectors. 

18. The combination of methods of claims 1. 9 and 17 for supporting the Ritz- 
vectors-based inverse numerical acoustics engineering process in 
particular but not limited to automotive engine and powertrain noise 
radiation application. 

19. The combination of methods of claims 1 and 9 for supporting the 
evaluation of multiple design alternatives, with respect to the acoustic 
performance prediction of vibrating structures. 

20. The combination of methods of claims 1. 9 and 12 for supporting the 
evaluation of multiple design alternatives, with respect to the modal-based 
acoustic performance prediction of vibrating structures. 

21. The combination of methods of claims 1 and 9 for supporting an 
automated structural-acoustic optimization process, with the objective 
function being the acoustic performance prediction of vibrating structures. 

22 'I^o^I^ ina ii 0n *. 0f , meth0d . 3 ° f C,a,ms 1 - 9 and 12 for supporting an 
structural-acoustic optimization process, with the objective 

!EKiiS2,r£ e m0da, - based acoustic Performance prediction of 



29 32 IB 480528 p. 26 R-592 Job-616 

13: 00 BIRD & GOEN & CO 32 16 480528 P. 26/35 



23 



23. ^processing engine adapted to carry out any of the methods of claim 1 to 
24 onS , ^ Uter w Pr0 ? ram prDdUGt for execut,n 9 on a computer, the computer 

Zm^SSSS any of the method steps of c,aim 1 to 22 

25 Vo^ th ? d ° f i nputtin 9 at a near terminal a representation of a body and 
coord.nates of a reference point and transmitting these to a remote 

Sr^JS r T Wn9 * f r ° 9ram for executi "9 any of the methods of claim 1 o 
22. and receiving at a near location an output of any of the methods. 

26. The method according to claim 25, wherein the output is one of: 

surfac^Vibra^n 'SfS? T^"** SUCh * S an ac ° ustic P ressure a 
body V,brat,on of bod V. a rev 'sed design of at least a part of the 



03-07-00 11:29 32 16 480528 P. 05 R-592 Job-616 

03-07-2000 12=48 BIRD & GOEN & CO 32 16 480528 P. 05x35 



ABSTRACT 

The invention relates to Computer-Aided Engineering (CAE) systems It 
concerns a new methodology to predict (1) the acoustic radiation 
charactenstics of a mechanical structure, under operational conditions, and 
(2) to identify the sources on a vibrating structure from measured sound 
pressure levels in the field. The methodology is based on a new approach to 
evaluate acoustic transfer vectors, based on the reciprocity principle and 
combined with interpolation techniques. 



The same methods are applicable to other vibrating energy forms which can 
be descnbed by the wave equation such as electromagnetic waves. 
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Figure 1: Acoustic Radiation Prediction. 
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Field Sound Pressure 




Figure 2: Inverse Numerical Acoustics/ 
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p M = <ATV(co)> {v n (g>) } 




Figure 3: The ATV-concept 
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Figure 4: Direct method for evaluating ATVs. This has to be repeated n 
times. 




Figure 5: Reciprocal method for evaluating ATV's. 
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100- ISO. 200. 250. 300. 350. 400. 



Figure 6: linear interpolation (in blue) of an Acoustic Transfer Vector (in red). 
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Figure 7: Spline interpolation (in blue) of an Acoustic Transfer Vector (in red). 




Figure 8: Panel Acoustic Contribution Analysis using a MDDBEM model 
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Figure 9; Optimization process based on ATVs. 
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